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OPTIMAL  DESIGN  OF  MULTIPLE  TUBE 
GRAIN  AND  THEIR  ARRANGEMENT* 


Zhao  Baihua 


ABSTRACT 


This  paper  presents  four  arrangement  methods  of 
multiple  tube  grain  in  the  combustion  chamber  and  the 
theorectic  calculating  formulae  for  the  grain  number 
n  and  relative  exterior  diameter  D.  Several  formulae 
for  computing  grain  number  n  in  the  "hexagonal  crowded 
arrangement"  have  been  corrected.  Results  of  the  com¬ 
puter  calculations  have  been  optimized.  Finally,  a 
group  of  grain  design  curves  have  been  given  to  meet 
the  requirements  of  various  web  coefficient  e.  These 
curves  have  practical  significance  to  the  grain  de¬ 
sign  of  some  missiles,  rocket  launchers  and  rocket 
boosters. 


*  Received  on  September  12,  1980 
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W. 


sp 


P 

'Yp 

Di 


T__d_ 
d"  Dj 


6"  Di 


e 


exterior  diameter,  interior  diameter,  length 
and  thickness  of  grain i 

grain  number j 

number  of  circle  arranged) 

total  impulse) 

specific  impulse) 

weigit  of  grain) 

specific  weight  of  propellant) 

interior  diameter  of  combustion  chamber) 

relative  exterior  diameter  of  grain) 

relative  interior  diameter  of  grain) 

web  coefficient) 

inlet  parameter) 
interior  inlet  parameter) 
exterior  inlet  parameter) 

ratio  of  interior  inlet  parameter  to  exterior 
inlet  parameter) 

packing  coefficient 
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FOREWORD 


Certain  missiles,  rocket  launchers  and  rocket  boosters 
usually  require  instant  completion  of  combustion  of  multiple 
tube  grain  in  order  to  achieve  high  initial  velocity.  However 
the  grain  number  n  in  currently  available  charts  and  tables  is 
less  than  40,  and  is  limited  to  only  one  arrangement  -  hexa¬ 
gonal  crowded  arrangement.  Furthermore  there  are  some  doubts 
as  to  some  of  the  theoretic  formulae  in  computing  the  value 
of  n  in  such  hexagonal  crowded  arrangement.  Thus  it  is  rather 
difficult  in  quickly  arriving  at  a  set  of  design  when  the 
grain  number  is  more  than  50. 

This  paper  performs  research  on  the  arrangement  methods 
of  multiple  tube  grain  in  combustion  chamber,  suggests  four 
feasible  arrangements,  and  derives  formulae  in  computing  n 
and  D  in  several  arrangements,  thus  correcting  several  for¬ 
mulae  used  in  computing  n  in  hexagonal  crowded  arrangement. 
Based  on  the  control  formulae  of  n  and  D  in  various  arrange¬ 
ments  and  the  control  formulae  that  relate  design  parameters 
and  geometric  parameters,  computer  programs  are  written  and 
large  quantity  of  calculation  results  are  obtained.  Finally 
through  selection  a  set  of  design  curves  having  a  range  of 
n  from  3  through  200  and  meeting  the  requirements  of  various 
e,  are  compiled.  Optimal  design  can  be  quickly  arrived  at  by 
using  the  methods  and  design  curves  presented  in  this  paper. 


METHODS  OF  ARRANGEMENT 


There  are  primarily  four  arrangement  methods  of  multiple 
tube  grain  in  an  engine*  namely  hexagonal  arrangement*  fixed 
peripheral  arrangement,  revised  arrangement  and  even  arrange¬ 
ment.  The  following  is  a  brief  discussion. 

1.  Hexagonal  Arrangement  (represented  by  ©) 

Hexagonal  arrangement,  otherwise  known  as  hexagonal  crowded 
arrangement,  is  an  arrangement  in  which  one  tube  is  first  placed 
in  the  center  of  a  combustion  chamber  and  then  other  tubes  are 
placed  hexagonally  outward.  Only  one  l/6  sector  is  being  studied 
here  in  Fig.  1  as  the  six  sectors  are  symetrically  identical. 
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By  progression  analysis*  when  the  number  of  circles  arranged 
i  is  positive  integer,  the  computing  formula  for  n  is 

n  ®  { 1  +3  I  (  1  +  I  »1 

Analysis  indicates  that  n  in  hexagonal  arrangement  depends 
on  both  the  value  of  i  and  the  odd/even  quality  of  i.  Based  on 
the  concept  of  sets*  the  computing  formula  for  n  in  hexagonal 
arrangement  can  be  represented  in  series  as 

*i,(  i  )  {a„  {a,  4- 8  k  }*.»,>{  O) 

CU=  1  +3  i  (1  +  «  ) 

.  (  { k  I  k  “2m—  l  <  i }  when  i  is  odd  number 

Combining  «,=■! 

(  i }  when  i  is  even  number 

m  is  natural  number  1,  2,  3,  4,  5,  6,  7  . 

All  values  of  n  in  hexagonal  arrangement  can  be  accurately 
computed  by  using  Formula  (1),  and  the  results  agree  with  the 
analysis  from  Fig.  1,  as  given  in  Table  1  where  dn  is  the  di¬ 
fferential  in  grain  number. 

Table  1. 


When  i  is  a  positive  integer*  the  relative  exterior  diame¬ 
ter  D  can  be  computed  by  using  the  following  formula 

°m 2  7 +  1 

Other  values  of  D  can  be  obtained  by  geometric  relation. 


I  came  across  several  formulae  In  computing  n  in  hexagonal 
arrangement  s 


1  +  3  (m  +■  m‘) 

(2) 

1  +  3  C(  m  —  1  )  +  (m  “  3  )  +  (»*•  —  5) 

+  (m-  7  )  +  •••} 

(3) 

1  +6m 

(4) 

where  m  is  natural  number.  There  are  some  doubts  about  these 
formulae.  Formula  (2)  is  presented  in  tl)  in  the  bibliography. 

The  n  values  computed  from  this  formula  are  by  majority  medium 
packing  density,  and  all  values  of  n  of  optimal  high  packing  den¬ 
sity  are  dropped.  Similarly  Formula  (3)  drops  a  good  part  of 
optimal  n  values j  some  nonexistent  n  values  are  present.  Formula 
(4)  simply  treats  hexagonal  arrangement  as  a  variation  of  arith¬ 
metical  progression  of  multiple  of  six,  and  it  is  obviously  an 
error.  Formula  (1)  suggested  in  this  paper  can  correctly  compute 
all  n  values  in  hexagonal  arrangement. 


2.  Fixed  Peripheral  Arrangement  (represented  by  A) 

Fixed  peripheral  arrangement,  otherwise  known  as  fixed-peri¬ 
pheral-constant-diameter-inward  arrangement,  is  an  arrangement  in 
which  integral  number  of  grain  sure  first  placed  on  the  most  outer 
ring  inside  a  combustion  chamber,  and  then  other  tubes  of  equal 
diameter  are  placed  in  circles,  inwardly  toward  the  center.  Ana¬ 
lysis  reveals  that,  for  the  same  values  of  n,  in  some  cases  the 
fixed  peripheral  arrangement  results  in  higher  packing  density 
than  the  hexagonal  arrangement.  For  instances 

Let  us  assume  n^&10)  being  the  grain  number  on  the  roost 
outer  ring,  the  number  of  eccentric  circles  being  i,  m  being  a 
natural  number,  we  can  obtain  the  following  formula  for  n  based 
on  the  geometric  relation  in  Fig.  2. 
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i  )=niT  2j" 

o  ! 

i 

•  =  i 

arc  sin  —  (2  *  -1  )D  l.A 

i 

j  (5) 

i 

n 

8inlT 

1  +sm— 

"i 

J 

When  n2*1^5»  no=0» 

when  n2*6~9»  no=l» 

when  n2>9.  then  nQ«0» 

n2  being  the  most  inner  circle  and  zh  being  the  sum  based  on  posi¬ 
tive  integers.  Formula  (5)  can  be  used  to  compute  the  grain  number 
on  each  eccentric  circle  as  well  as  the  total  grain  number  in  the 
fixed  peripheral  arrangement. 


Fig.  2.  Diagram  of  fixed  peripheral  arrangement. 


3.  Revised  Arrangement  (represented  by  X) 

Revised  arrangement  is  an  arrangement  in  which  integral  num¬ 
ber  of  grains  are  first  placed  in  the  most  inner  circle  inside  a 
combustion  chamber  (but  this  is  not  hexagonal  arrangement),  and 
then  other  tubes  of  equal  diameter  are  placed  in  circles,  toward 
the  side  of  the  chamber.  It  is  an  arrangement  derived  from  revi¬ 
sing  the  fixed  peripheral  arrangement  and  can  in  some  cases  achieve 
higher  packing  density  than  the  fixed  peripheral  arrangement. 

Let  us  assume  n2(>3)  being  the  grain  number  on  the  most  inner 
circle,  the  number  of  eccentric  circles  being  i,  we  can  arrive  at 
the  following  formula  for  n_(i)  based  on  the  geometric  relation  in 
Pig.  3. 
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«.( » )-»,+  s 


arc  sm 


1  -Di*i  U 


1  +(2  »'  —  1  )  sin- 


When  n2*3,4,5,  nQ*Os  when  ng *6, 7,8,9,  nQ*l}  zh  is  the  sui, 
based  on  positive  integers.  Di_  is  the  relative  exterior  dia- 
meter  of  the  ith  circle  when  the  grain  number  on  the  most  inner 
circle  is  n„ 


Fig.  3.  Diagram  of  revised  arrangement. 


Formula  (6)  can  be  used  to  compute  the  grain  number  on  each 
eccentric  circle  as  well  as  the  total  grain  number  nx ( i )  in  the 
revised  arrangement. 


4.  Even  Arrangement 


Experiments  indicate  that  for  the  above  arrangements,  even 

*  I 

when  relatively  high  dispersion  of  internal  ballistic 

performance  will  still  exist,  causing  dispersion  of  initial  ve¬ 
locity  in  rockets  and  missiles.  Such  phenomenon  in  crowded  arran 
gement  is  caused  by  fracture  of  grain  tubes  due  to  partial  flow 


velocity  differential  and  pressure  differential  created  m  the 
interior  and  exterior  inlet.  In  order  to  solve  the  problem  in 
dispersion  of  internal  ballistic  performance  and  grain  tube  sta¬ 
bility,  such  even  arrangement  is  presented. 

_ln  Dp-*Dth  (7) 

Dth  the  theoretic  value  of  the  relative  exterior  diameter 
of  the  grain,  Dp  is  the  practical  value  of  the  relative  exterior 
diameter  of  the  grain  for  revision  coefficient  $  for  consideration 
As  far  as  consideration  limiting  case  for  crowded  arrangement  is 
concerned,  <j>*0.98  is  usually  adopted!  however  for  even  arrangement 
<j>=0.95,  0.90,  0.85  are  employed.  Table  2  lists  the  value  of  n  and 
the  computed  results  of  its  D  for  several  arrangements. 


Table  2.  D+.  *f.  (n,  arrangement) 


arr/  Dth 


X 

0.1137 

® 

0.1213 

A 

0.1155 

X 

0.1149 

A 

0.1114 

® 

0.1111 

X 

0.1039 

A 

0.1076 

X 

0.1074 

A 

0.1040 

® 

01029 

A 

0.1007 

X 

0.0985 

A 

0.0976 

X 

0.0960 

108  !  ® 


181  ®  j 

123  I  x  I 
12S  |  A  : 
127  |  © 


139  '  ® 


Two  sets  of  control  formulae  are  combined  as  follows  in  order 
to  facilitate  computer  programming. 

n-fj^D,  arrangement,  i)  control  formulae 

»e(  »)-{«*>  {a,  +  6  D--2  —  - 

(The  other  formulae  are  omitted  here) 


«.(  «  )-« o+  2 


arc  sin  —■ — —  I 

1  —  Dim,  .* 


D  = 


It 


1  +(2  i  —  1  )sin  — 


Design  parameter=f2  (geometric  parameter)  control  formulae 

5*  ( “ »  A »  arrangement ) 

ti  =  «  (o2-<r)  =  n  ( » ,  X  *  arrangement) 

1^fe-=~iT(5_a)c  1  "  *  (5*-S*)D-Q(«.  X,  arrangement) 

«,  -  =«,(  n ,  x ,  arrangement) 

®p*^°th  (n’  arrangement),  $*0.98,  0.95.  O.90,  0.85 


Computation  results  are  compiled  in  the  table  that  follows. 


9 


OPTIMAL  DESIGN  OP  MULTIPLE  TUBE  GRAIN  AND  THEIR  ARRANGEMENT 


For  n  ranging  from  3  to  200,  there  are  21  values  of  n0(i)  in 
hexagonal  arrangement}  44  values  of  n^i)  in  fixed  peripheral  arr¬ 
angement  and  30  in  revised  arrangement.  To  facilitate  applications 

and  optimal  design,  five  categories  are  listed  to  the  size  of  their  packing 
coefficient  n  (<}>  +  0-98,  X  +  1).  The  classification  and  screening  results 
ate  shown  in  Table  3. 


Table  3* 


1  category 

arr. 

jn  values  (from  left  to  right  on  decreasing 

00.70) 

® 

109,  107,  151,  OS,  121,  163,  55,  109,  139 

A 

19 

1* 

<0.04—0.60) 

! 

© 

31,  7,  73 

mm 

173,  101,  132.  100,  130 

a 

101,  130,  130,  93,  90,  180,  190 

"E 

<0.00—0.04) 

© 

109,  19.  01,  91.  127,  97 

a 

119,  12S,  195,  ISO,  144,  165,  60,  101.  90,  30.  112,  95,  70,  106,  05,  41, 

74.  52,  50,  31,  20,  45,  34 

X 

153,  ISO.  123,  141,  60,  193,  17,  149,  07,  70,  102,  02.  57,  50,  40,  62, 

S3,  00 

1 

’ll  1 

1 

<0.5—0.60)  { 

© 

43 

A 

8,  44*  22*  l$o  13,  27,  24 

X 

J  47,  6,  3 

-  i 

«0.5)  1 

mm 

9,  10 

mm 

4,  5.  43 

The  relations  of  7)j ,  with  respect  to  n,  arrangement  and  \ 
are  obtained  through  selection  and  are  illustrated  in  Fig.  4  and 
Pig.  5. 

Analysis  through  categorized  selection  reveals  that  within  the 
range  of  there  are  9  arrangements  in  hexagonal  arrangement}  i.  e., 
n0*199»  187#  151#  85#  121,  163#  55#  109#  139  whereas  in  fixed  peri¬ 
pheral  arrangement#  there  is  only  one,  namely  n^»19.  Within  the 
range  of  *)x  there  are  3  arrangements  in  hexagonal  arrangement}  i.  e., 
n@  *31#  7#  73  and  five  arrangements  in  fixed  peripheral  arrangement} 
i.  e.#  nA«173#  181#  132#  188,  138  as  opposed  to  seven  arrangements  in 
revised  arrangement}  i.  e.,  nx*l6l,  136,  130,  93,  98,  180,  196. 
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ig.  4.  iQr“  *r)I(nt  arrangements)  pig.  5.  =  Oj.(n»  arrangements) 


In  order  to  satisfy  the  need  in  various  applications*  values 
of  n  are  divided  into  four  groups  and  are  further  selected  by  cate¬ 
gory.  Results  are  given  in  Table  4. 

Table  4. 


n  range 


Group  1 

8— SO 


Group  2 

51—100 


Group  3 


101—150 


Group  4 


ii 


hi 


’ll 


’ll 


*ii 


11 


n  values  (from  left  to  right  on  ^  ) 


19A 


31®,  7® 


19®,  38  A.  37®,  41  A.  50  X,  40  x ,  31  A,  20A.  I8A,  3«A 


85®,  55® 


73®,  93  x,  98  x 


50  A,  66  x  ,  90r;.,  87  x,  95A,  70A,  Tlx,  82  x ,  61®.  Ii3_,  57  x  ,  91®,  74A. 
12!®,  109®,  139® 


’ll 


136  x,  130  x,  I32A,  138A 


’ll 

’ll 


11®,  125A,  144  A,  123  x.  Mix,  I01A.  117  .  1 10  x  ,  112.'.,  106A,  102  x, 127® 


199®,  187®,  lol®,  183® 


151—200 


The  relationship  of  packing  coefficient  with  respect  to  n, 
arrangement  and  \  in  each  group  are  given  from  Fig.  6  through  Fig. 

9.  As  far  as  Ox  is  concerned,  there  is  one  in  group  1,  i.  e. ,  19A  i 
two  in  group  2,  i.  e.,  850,  55© 1  three  in  group  3»  i*  e. ,  1210,  109® 
1390,  and  four  in  group  4  i.  e.,  199®,  1870,  1510*  and  I630.  As  far 
as  Oi^s  concerned,  there  are  two  in  group  1,  i.  e.,  31©»  7©*  three 
in  group  2,  i.  e.,  73© •  93X,  98X1  four  in  group  3»  i.  e.,  136X,130X, 
132a,  138a,  and  five  in  group  4,  i.  e.,  161X,  173^-,  181a,  180X  and 
196X. 

Due  to  different  arrangement,  for  equal  values  of  grain  number 
n,  their  values  in  0  are  different  and  in  some  cases  the  difference 
is  quite  big.  For  example,  at  n*19  when  hexagonal  arrangement  is 
changed  to  fixed  peripheral  arrangement,  its  value  in  0  can  raise 
from  0jLto  "Ox  •  At  n*85»  OjiCan  also  be  raised  to  TQj  if  fixed  peri¬ 
pheral  arrangement  is  changed  to  hexagonal  arrangement.  Therefore 
we  can  say  that  an  optimal  arrangement  exists  when  values  of  n  are 
equal.  Optimal  design  for  n*13,  19,  31,  43,  85,  151  is  listed  in 
Table  5. 

Table  5*  Optimal  arrangement  for  same  values  of  n. 
value  of  n  arr.  n  optimal  arrangement 


The  above  analysis  shows  how  to  achieve  optimal  design  of  mul¬ 
tiple  tube  grain  and  their  arrangement  according  to  the  principle 
of  maximum  packing  density.  Certain  launchers  and  uncontrolled  roc¬ 
ket  boosters  demand  extreme  short  working  time  and  the  requirement 
is  very  stringent.  In  other  words,  there  is  a  definite  requirement 
on  the  value  of  e^.  The  following  analysis  will  present  optimal  de¬ 
sign  of  multiple  tube  grain  and  their  arrangement  in  the  cases  that 
the  values  of  ej  is  to  be  satisfied. 

When  e^O.l,  we  can  select  corresponding  grain  type  to  fulfill 
the  requirement  of  e^  of  various  ranges.  For  ¥^<0.1,  available 
grain  types  are  multiperforated  grain,  laminated  grain,  coiled 
grain  and  multiple  tube  grain,  the  first  being  an  area  augmenting 
grain.  Generally  multiple  tube  grain  is  used  in  order  to  achieve 
near  constant  area  behavior.  Analysis  shows  that  for  a  given 
there  are  multiple  values  of  n  and  its  arrangement  that  can  meet 
the  requirement.  Results  are  shown  in  Fig.  10  when  optimal  selec¬ 
tion  is  performed  based  on  maximum  packing  coefficient.  Analysis 
on  Fig.  10  reveals  the  foil owing i 

¥^*0.02«'0.044  is  the  region  of  high  packing  density  web  coe¬ 
fficient  e^.  In; such  e^  region  values  of  n  and  its  arrangement  in 
ranges  of  or  TQj  can  be  selected  which  will  satisfy  e^  provided 
that  pressure  is  controlled  properly,  “0.044*0.07  is  the  region 
of  medium  packing  density  web  coefficient  'e^.  In  such  region  the 
majority  of  n  values  selected  on  e^  fall  on  grange,  except  for 
a  few  cases  of  n  and  its  arrangement  based  on  n,  .  ¥^*0.07*0.11  is 
the  region  of  low  packing  density  web  coefficient  e^  In  such  re¬ 
gion  with  the  two  exceptions  of  e^^  (e^O. 708(1^)  and  e1«0.66l(ir)i)J 
a  great  majority  of  n  values  selected  on  e^  and  its  arrangement 
are  low,  falling  on  grange  or  range.  For  a  required  e^  whose 
corresponding  n  and  V)  from  its  arrangement  are  on  the  low  side, 
e^  value  can  be  modified  through  the  selection  of  a  propellant  of 
constant  combustion  rate  or  through  the  selection  of  pressure , 
thus  boosting  its  corresponding  n  and  *)  on  its  arrangement. 
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Fig.  10.  Relation  of>r)*f(e1)  when  \*1. 

Thus  we  are  able  to  select  an  optimal  n  and  its  arrangement 
to  meet  the  requirement  of  grain  packing  design. 

This  paper  constructs  the  following  three  sets  of  multiple 
tube  grain  design  curves  based  on  computation  results  with  res¬ 
pect  to  n  and  its  arrangement. 

d*d(n,  arrangement ,  \ ,  <J> )  design  curvet 
e,«e,(n,  arrangement,  7\,  <J>)  design  curve j 

w 

*Q(n,  arrangement,  design  curve. 

These  sets  of  design  curves  are  presented  in  Fig.  11  through 
Fig.  16.  By  employing  these  sets  of  multiple  tube  grain  design 
curves  an  optimal  design  of  multiple  tube  grain  and  its  arrangement 
satisfying  e^  requirement  can  be  quickly  obtained. 
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APPLICATION  EXAMPLE 


Suppose  that  a  certain  engine  has  an  interior  diameter  (S^> 
of  90  mm  and  a  total  impulse  (I+20pc^  of  *90  kg-sec  is  to  be  gua¬ 
ranteed  for  initial  velocity.  Try  to  design  an  optimal  tube  grain 
and  its  arrangement  for  a  combustion  time  t^ 0.018  sec  at  +20°C. 

[Solution]  Select  a  solid  propellant  whose  specific  impulse, 
specific  weight  and  combustion  rate  are  as  follows « 

ISp*242  second,  Yp*1.704  g/cnr*’  *6.15p°*^ 

Pressure  at  ♦20°C  in  combustion  chamber  p+2Oo0*3OOkg/cm 
Then  the  grain  thickness  is 

e1=2apntb«2X6. 1 5X300° *36XO. 018-1. 73(mm) 

The  required  web  coefficient  is 

® l=el^Dl*1 *  73/90-0.0192 

The  required  grain  weight  is 

I  190 

W  -242*0-785<k«> 

e^  for  consideration  is  very  small,  and  for  the  purpose  of 
minimizing  performance  dispersion  and  stabilizing  multiple  grain, 
even  arrangement  for  multiple  tube  grain  is  chosen.  Select  $-0.90, 
for  e^ -0.0192,  from  Fig.  15,  at  $-0.90,  the  following  ten  propo¬ 
sitions  are  obtained  as  listed  in  Table  6  through  the  design 
curve  e-e^n,  arrangement,  K). 


Table  6. 


i 

55® 

«SA 

78® 

82  X 

15® 

S3® 

m© 

0.80 

1.08 

1.025 

OSS 

1.15 

1.175 

1.085 

1.01 

r)  category 

I 

n 

D 

E 

8 

I 

m 

m 

I 

I 

First  look  at  the  four  propositions  belonging  to  7^,  i.  e., 
n-55©*  85©,  1090  and  1210.  From  Fig.  4  it  can  be  seen  that 

T)n«85  y)n-121  >*)n«55  V)n-i09*  looking  up  Fig.  16  for  those  four 
cases  belonging  to *]jj,  f«0.90,  Wp/y>pD?k*ft(n,  arrangement,  \) , 

18 


four  values  of  n  and  \,  il,  *  corresponding  to  its  arrangement 
are  obtained,  as  indicated  in  Table  7* 


Table  7. 


II 

5S© 

35® 

1 

100® 

131© 

X 

0.M 

0.09 

1.935 

1.05 

Q 

0.0047 

0.0040 

0.00335 

0.0037 

X 

131 

150 

104 

173 

Since  such  propellant  has  a  threshold  ***150  and  the  pressure 
for  consideration  are  relatively  high  (300kg/cm  ),  the  allowed  x* 
has  to  be  raised  accordingly,  and  thus  the  proposition  of  n*85 
hexagonal  arrangement  is  selected. 

At  n=85©  and  \=0.89,  from  Fig. 14  we  get  d=0.050.  From  Table 
2  we  get  D=<J)Dth=0. 9X0. 0984*0. 08856,  Finally  we  obtain  the  optimal 
design  that  meets  the  requirement  of  armament  technology  asi 

— LXn  *  -  136  X  85  (hexagonal  arrangement) 

4  4.50 

After  rounding  we  get 

-2- -LXn  (arrangement)  *  137  X  85  (hexagonal  arrangement). 

d  4.6 
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